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We estimate the consequences of finite masses of pseudoscalar mesons on the decay rates of 
scalar glueballs in the Witten-Sakai-Sugimoto model, a top-down holographic model of low-energy 
QCD, by extrapolating from the calculable vertex of glueball fields and the p' meson which follows 
from the Witten-Veneziano mechanism for giving mass to the latter. Evaluating the effect on the 
recently calculated decay rates of glueballs in the Witten-Sakai-Sugimoto model, we find a strong 
enhancement of the decay of scalar glueballs into kaons and rj mesons, in fairly close agreement with 
experimental data on the glueball candidate /o(1710). 
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The fundamental theory of the strong interactions, 
quantum chromodynamics (QCD), which has quarks con¬ 
fined in color-neutral bound states, admits also bound 
states whose valence constituents are all gluons, the non- 
abelian gauge bosons of QCD. This prediction of addi¬ 
tional mesons called gluonia or glueballs dates back to 
the eaidy 1970’s [l| and has been substantiated by lattice 
QCD y , which estimates the mass of the lowest glueball 
state to be around 1600-1800 MeV. Experimentally, how¬ 
ever, their status remains unclear and controversial Q. 
The lowest scalar glueball state has quantum numbers 
of the vacuum and can be expected to mix with scalar 
mesons made from quarks and antiquarks. To disentan¬ 
gle the contributions, information on decay processes is 
needed. Theoretical expectations vary greatly—the low¬ 
est glueball state may be even so broad that it forms a 
mere background for the isoscalar meson spectrum [^. 

QCD in the limit of a large number of colors A^c [EQ, 
which in many cases turns out to be a remarkably suc¬ 
cessful approximation to real QCD with = 3, predicts 
a parametric suppression of decay rates of glueballs com¬ 
pared to light quarkonia by a factor 1/iVc as well as a 
suppression of mixing Q. If glueballs are indeed narrow 
and not strongly mixed, one should be able to identify one 
of the isoscalar-scalar mesons below 2 GeV as a predom¬ 
inantly glueball state. In phenomenological studies the 
experimentally well-established Q mesons /o(1500) and 
/o(1710) have been identified alternatingly as possible 
glueball candidates [§-[ll|. Both are comparatively nar¬ 
row states, but their decay patterns are rather different: 
/o(1500) decays primarily into four pions and secondly 
into two pions, with decays into kaons and rj mesons 
suppressed, whereas /o(1710) instead decays predomi¬ 
nantly into two kaons, with a ratio @ T(2t:) /V{KK) = 
0.41 ([ij, much lower than 3:4 expected from a flavor- 

blind glueball. In the case of /o(1500) the strong devi¬ 
ation from flavor-blindness is usually attributed to mix¬ 
ing, while for /o(1710) it has been suggested that glue¬ 
balls couple more strongly to the more massive pseu¬ 
doscalar mesons, a mechanism termed “chiral suppres¬ 


sion” 12, 1^, which could make it possible that /o(1710) 
is a nearly unmixed glueball as most recently argued for 

in [13,113 (see also Ref. 0). 


Since lattice QCD results on glueballs in interaction 
with quarks are still sparse, in particular concerning de¬ 
cay patterns, it is of interest to employ (top-down) gauge- 
gravity duality, a string-theoretic approach to study 
strongly coupled large-fVc gauge theories, to obtain new 
insights from first principles 15|. In fact, the spectrum 
of glueballs has been one of the first applications of a 
nonsupersymmetric holographic model derived by Wit¬ 
ten [l6| from type-IIA superstring theory 17H19j. The 
Witten model has subsequently been extended by Sakai 
und Sugimoto to include chiral quarks through D8-D8 
probe branes [i^, IHj . With only one free coupling con¬ 
stant at a (Kaluza-Klein) mass scale Mkk ~ 1 GeV, this 
provides a remarkably successful model for low-energy 
QCD with quantitative predictions for vector and axial 
vector meson spectra and decay rates that agree with real 
QCD to within 10-30% 

In Ref. [2^, the Witten-Sakai-Sugimoto model was 
used for the hrst time to evaluate the decay rate of the 
lowest glueball state into pions and to compare with ex¬ 
perimental data for the /o(1500), although the mass of 
the lowest hol(OTaphic glueball mode comes out at 855 
MeV. In Ref. 2^ we have recently revisited this calcu¬ 
lation with the result that the decay width of the lowest 
mode is much higher than the one obtained in . Since 
the lowest mode corresponds to an “exotic polarization” 
[l3| of the gravitational fields of the Witten model, we 
have proposed to discard the latter and to instead con¬ 
sider the next-lowest, predominantly dilatonic mode with 
mass 1487 MeV as corresponding to the glueball in QCD. 
Despite the closeness of its mass to that of the /o(1500) 
meson, we have found that the decay pattern into two 
and four pseudoscalar mesons is not reproduced: the de¬ 
cay rate of /o(1500) into two pions is underestimated by 
about a factor of 2, while the prediction for the domi¬ 
nant decay mode of /o(1500), which is decay into four 
pions, is an order of magnitude too small. Extrapolating 
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the mass of the holographic glueball to that of the glue- 
ball candidate /o(1710) (which is within 16% of the mass 
of the dilatonic mode) [^, we have instead found close 
agreement with the decay rate into two pions. Since the 
Witten-Sakai-Sugimoto model is chiral, pions, kaons, and 
rj mesons are predicted simply in flavor-symmetric ratios 
3:4:1, thus failing to explain the much stronger decays 
into two kaons and two rj mesons. 

In this Letter we study the possible effect of finite 
quark masses on the decay rate of glueballs in the Witten- 
Sakai-Sugimoto model in order to see whether a suffi¬ 
ciently strong enhancement of the decay into kaons and 
also rj mesons could result. In Refs. [n, m , it has 
been shown that nonlocal mass terms implementing Gell- 
Mann-Oakes-Renner relations can be induced by either 
worldsheet instantons or a deformation by a bifundamen¬ 
tal field related to the open string tachyon that arises 
between (parallel) D and anti-D branes. Since no com¬ 
plete calculation along these lines exists, the additional 
coupling of glueballs and pseudoscalar mesons induced by 
the nonlocal mass term is not known. However, the addi¬ 
tional coupling of glueballs to rj' mesons due to the part of 
its mass term which arises from the anomalous breaking 
of U(1 )a flavor symmetry can be calculated exactly. We 
propose to use that as a simple and also plausible model 
for how the totality of nonlocal mass terms for the pseu¬ 
doscalar meson depend on the scalar glueball fields, and 
to thereby extrapolate the results for the glueball decay 
pattern obtained in [2^ to finite pseudoscalar masses. 

In the chiral Witten-Sakai-Sugimoto model 0, El, 
the coupling of pseudoscalar Goldstone bosons as well 
as vector and axial vector mesons to scalar and tensor 
glueball fields are determined by the dependence of the 
Dirac-Born-Infeld part of the D8 brane action on metric 
and dilaton fields. 


= -TosTr J v"- det {~gMN + 2Wi^Mw), 

( 1 ) 

where $ and qmn are the dilaton and the nine¬ 
dimensional induced metric in the ten-dimensional back¬ 
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with circle-compactified x^ x'^ -\- 27r/MKK and Mkk = 
The stacks of Nj D8 and anti-D8 branes 
are assumed to be localized at antipodal points, giving 
rise to trivial embeddings x^ = const., which extend from 
the holographic boundary at u = oo to the minimal point 
ukk where branes and antibranes connect. This breaks 
the chiral group U(iV/)j;, x U(7 V/)h down to its diagonal 


group, leading to a nonet (for Nf = 3) of pseudoscalar 
Goldstone bosons described by 

/ OO 

dzH^(z,a:) = (4) 

-CX3 


where zjuKK = \/(u/ukk)^ — 1 parametrizes the radial 
extent of the joined D8 and anti-D8 branes. Matching 
the pion decay constant to 92.4 MeV and the mass of 
the lowest vector meson mode Af^{z,x) = p^{x)ip(i){z) 
to the p meson mass nip ~ 776 MeV fixes Mkk = 949 
MeV and A = gyyi^c = 16.63 [2^, El; matching in¬ 
stead THplsja with a the st ring tension of the model to 
large-Vc lattice results [H, l24| gives a somewhat lower 
value of the ’t Hooft coupling, A = 12.55, which we use 
with the higher value to give a band of variation for the 
holographic predictions. 

In Ref. [ 2 ^ (see also [ 2 ^) it was shown that 
the U(I)^ anomaly requires to combine the Ramond- 
Ramond 2-form field strength F-i with the isoscalar tjq 
that is localized on the D8 branes in a gauge-invariant 
combination F 2 with bulk action 

with 


F 2 


67ru|^K-^KK 





du A dx '^, 


( 6 ) 


where 6 is the QCD theta angle and 


Vo{x) = j dzTTA^{z,x). 


(7) 


This gives rise to a Witten-Veneziano [E,E1 mass term 
for r]Q that is local with respect to the effective 3-1-1- 
dimensional boundary theory but nonlocal in the bulk, 
with mass squared 


ml = ^^A"M2 . 


( 8 ) 


For Nf = Nc = ^, Mkk = 949 MeV, and A varied from 
16.63 to 12.55 one finds toq = 967... 730 MeV. 

The other pseudoscalar mesons described by 0 are 
massless in the Witten-Sakai-Sugimoto model. Gurrent 
quark masses can in principle be introduced through a 
deformation by a bulk field T in the bifundamental rep¬ 
resentation of the chiral symmetry group that is re¬ 
lated to tachyon condensation or alternatively through 
worlsheet instantons [ 23 . Both introduce nonlocal mass 
terms for the pseudoscalar mesons, which one may qual¬ 
itatively write as 


f d'^x f duh{u)TT (r(u)Pe-*7<i^MG.x) \ ^ 

UKK 


( 9 ) 
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where h{u) includes metric fields. Choosing appropriate 
boundary conditions for T, the quark mass matrix arises 
through 


du h(u) T{u) (X A4 = diag(mu, md, mg), (10) 


thereby realizing a Gell-Mann-Oakes-Renner relation. 

Integration over u leads to mass terms for all Goldstone 
bosons including one for the flavor singlet 770 in addition 
to the Witten-Veneziano mass term. The flavor octet rjs 
and ? 7 o can be diagonalized to mass eigenstates rj and 
rj'. With Ai = diag(TO, m, mg), m = {iriu + md)/2, fixing 
= 140 MeV and rriK = 497 MeV, this diagonalization 
yields for A = 16.63 ... 12.55: 

= 518... 476 MeV, = 1077... 894 MeV, (11) 
6»p =-14.4°...-24.2°, (12) 


with 9p the octet-singlet mixing angle, which shows that 
the above holographic result for mp is in the right ball¬ 
park [^ . 

In order to determine how mass terms affect the cou¬ 
pling of glueballs and mesons worked out in Ref. [ 2 ^, 
we would need to know the dependence on dilaton and 
metric fields of h{u) as well as the profile of the bifun¬ 
damental field T{u). Absent those, we turn to the fully 
known nonlocal mass term produced by ([5]) and ([6]) for 
r]Q. Inserting the mode exp ansion of glueball fields Go 
and Ge defined in Ref. [2^, we find 


^vo' ^ J (1 - 3doGp -b ScpGs) -b ... 

(13) 

with 

17 915 

^ 33 _ 4 ^ 15.829 

“ A./^JV.Mkk ’ 

where the latter is given for completeness only, since we 
are going to discard the “exotic” mode Ge given the 
results in [ 2 ^. Here He^e{u) are the radial profile func¬ 
tions of the glueball modes, normalized such as to give a 
canonical kinetic term for Gd,e{x). 

Given the similarity to how a nonlocal mass term 
is generated through worldsheet instantons, this result 
seems to be a reasonable first guess of how nonlocal 
mass terms couple in general. As far as a bifundamen¬ 
tal field T associated with tachyon condensation is con¬ 
cerned, a plausible guess would be that the metric depen¬ 
dence derives from the integration measure of D8 branes, 
d^xe~^y/^. For the predominantly dilatonic glueball 
field [s^ , this turns out to have exactly the same depen¬ 
dence on terms linear in Gd{x)He{u) as follows from ([5]) 
and namely a factor {1 — 2,Gd[x)He(u)). In order to 
calculate the coupling constant analogous to dp in (fTTl) . 


one would need to know the holographic profile of T, of 
which we only know that it will be concentrated around 
u = mkk- As a simplistic guess one could try a func¬ 
tion that mimics the profile of the term in the 

D-brane action when A^ equals the zero mode describing 
the Goldstone bosons. This would simply determine the 
analog of dp to be equal to the coupling di appearing in 
the chiral Gdtttt term [ 2 ^ . 


where di ~ 17.226A“^/^V“^Mj^^. This differs from dp 
by a mere 4%. 

We shall therefore continue with the working hypoth¬ 
esis that the overall coupling of the glueball field to the 
mass term for the pseudoscalar mesons is universal. This 
essentially assumes that the mixing of singlet and octet 
mesons is invariant under a holographic renormalization 
group evolution, which in particular implies the absence 
of a direct coupling Gryry'. 

With this assumption, i.e. adding 

(17) 


with Pi the mass eigenstates of the pseudoscalar mesons 
to (fTB)) . we obtain the following modification factor for 
the decay rate to two pseudoscalar mesons of mass mp: 




M2 J 


1 -j- ct 


M2 


(18) 


with 


a = 4(3do/di - 1) = 8.480 ... for Go- (19) 


An analogous calculation for the “exotic” scalar glueball 
using the results of 2^ for the chiral contributions gives 


a = 4(5co - ci)/(ci -b 2 ci) = 2.630 ... for Ge- (20) 

[Note that to leading order the dependence on A and 
Mkk drops out in (HU) and (Uni).] In (113) the first fac¬ 
tor represents a simple kinematical suppression which is 
overcome by the coupling of the glueball field to the mass 
term of the pseudoscalar fields. A similar result, but with 
a = 1 was obtained in Ref. 3^ for a simple effective field 
theory where the scalar glueball field is identified with the 
dilaton of QCD (a scalar field with a potential matched to 
the QGD trace anomaly). With a = 1 the nonchiral en¬ 
hancement is cancelled by the kinematical suppression to 
order mp/M2, thus restoring approximate flavor symme¬ 
try, while for larger values mp the net effect is a (slight) 
reduction of the decay rate. 

In Table H] we compare the deviations from flavor sym¬ 
metry as they are reported by the Particle Data Group 
0 with the modification resulting from (1181) and (113. 
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Remarkably, the experimental ratio T{'K'k)/T{KK) is re¬ 
produced within the experimental error bar, whereas the 
prediction for T{'q'q)/T{KK) remains within 1.33 stan¬ 
dard deviations. 

In Table m we compare our complete set of predic¬ 
tions for the decay rates for a scalar glueball with mass 
corresponding to either 1505 MeV (/o(1500)) or 1722 
Mev (/o(1710)), with and without inclusion of masses 
for the pseudoscalar mesons, to experiment. In the case 
of /o(1500), the (experimentally well-known) decay pat¬ 
tern is neither matched qualitatively nor quantitatively 
0 . Inclusion of the pseudoscalar masses helps for the 
total width but modifies the decay pattern adversely. 
For /o(1710), branching ratios are less accurately known. 
The prediction for the total width is increased slightly 
above the experimental value when masses are included, 
but the decay pattern into two pseudoscalar mesons is 
improved markedly, as we have already shown in Table 
m In fact, the experimental results quoted for the partial 
widths should be considered as upper values as they ig¬ 
nore the possibility of decay into four or more pions, since 
at least decay of /o(1710) into two w mesons and further 
to six pions has been seen Q in J/ij) —>■ 7/o(1710) (in fact 
at the level of 75% of the rate into two pions 0, so that 
the holographic prediction may not be very far off) 

The only remaining major mismatch between existing ex¬ 
perimental data for /o(1710) and the prediction of 
the Witten-Sakai-Sugimoto model thus appears to be the 
rather high rate for decay into four pions, which is pre¬ 
dicted by the latter to proceed through 2 p and pirn at 
the level of about twice the rate into two w mesons 24| . 
and this prediction is not modified by the introduction 
of quark masses at the level of our approximation. 

To summarize, by extrapolating the exactly calculable 
coupling of scalar glueballs to the mass term of the isos¬ 
inglet pseudoscalar meson in the originally chiral Witten- 
Sakai-Sugimoto model, we found a significantly enhanced 
decay of scalar glueballs into kaons and rj mesons com¬ 
pared to flavor-symmetric ratios. This is in line with 
the previously proposed mechanism of “chiral suppres¬ 
sion” of scalar glueball decay which has been posited as 
explanation how the isoscalar meson /o(1710) with its 
preferred decay into two kaons could be predominantly 
gluonic rather than an ss state Elil. From this we 
conclude that the top-down holographic Witten-Sakai- 
Sugimoto model may well be consistent with a glueball 
interpretation of /o(1710) while disfavoring the other 
popular glueball candidate /o(1500). In this case, the 
successful reproduction of the branching ratios given in 
Table I is correlated to a sufficiently small rate for the 
decay G —>■ ryry', for which only upper limits exist so far 


37| . Moreover, the Witten-Sakai-Sugimoto model pre¬ 


dicts significant partial widths for the decay of /o(1710) 
into four and six pions, of which so far only the latter 
have been confirmed experimentally according to Q. 

We thank Denis Parganlija and Timm Wrase for many 


/o(1710) 

exp.(PDG) WSS massive 

|■^inn)/^iKK) 0.55 

0.463 

4:-r{pp)/r{KK) 

1.92±0.60 

1.12 


TABLE I. Flavor-asymmetric deviation of branching ratios 
of glueball candidate /o(1710) compared to the nonchiral 
enhancement in the Witten-Sakai-Sugimoto model resulting 
from (fTSll and (fT^ with mp = m,r,if,77 = {140, 497, 548} MeV. 
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compared to the prediction obtained in Ref. from mode Go in the chiral Witten-Sakai-Sugimoto model, and finally to the 
extrapolation of the massive case proposed in this Letter. 
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